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It is shown that the corrections of a 2 order to the total cross sections for interaction of elementary 
hydrogen- like atoms with target atoms, reported in the previously published paper [S.Mrowczyriski, 
Phys.Rev. D36, 1520 (1987)], do not include some terms of the same order of magnitude. That 
results in a significant contribution of these corrections in particular cases. The full a 2 -corrections 
have been derived and it is shown that they are really small and could be omitted for most practical 
applications. 
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I. INTRODUCTION 

The experiment DIRAC §, now under way at PS 
CERN's, aims to measure the lifetime of hydrogen-like 
elementary atoms (EA) consisting of tt + and n~ mesons 
(A2ir) with an accuracy of 10%. The interaction of tt + it~- 
atoms with matter is of great importance for the experi- 
ment as dissociation (ionization) in such interactions 
is exploited to observe and to measure its lifetime. In 
the experiment the ratio between the number of 7r + 7r~- 
pairs from A^ dissociation inside a target and the num- 
ber of produced atoms will be measured. The lifetime 
measurement is based on the comparison of this experi- 
mental value with its calculated dependence on the life- 
time. The accuracy of the cross sections for interaction 
of relativistic EA with ordinary atoms, which are behind 
all these calculations [|| , is essential for the extraction of 
the lifetime. 

Study of interactions of fast hydrogen-like atoms with 
atoms has a long history starting from Bethe. One of 
the resent calculations for hydrogen and one-electron ions 
was published in . Interactions of various relativistic 
EA consisting of e ± , tt ± , fi ± , were considered in 



to the EA total interaction cross sections obtained in 
(Through this paper a is the fine-structure constant/ 



II. GENERAL FORMULAS 

As shown in |?j analysis of the relativistic EA inter- 
action with the Coulomb field of target atoms can be 
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performed conveniently in the rest frame of the projec- 
tile EA (anti-lab frame). As the characteristic transfer 
momentum is of the order of the EA Bohr momentum, in 
this frame after the interaction EA has a non-relativistic 
velocity and thus the initial and final states of EA can be 
treated in terms of non-relativistic quantum mechanics. 
In this manner the well-known difficulties of the relativis- 
tic treatment of bound states can be get round. 

As in the EA rest frame a target atom moves with the 
relativistic velocity, its electromagnetic field is no longer 
pure Coulomb. It is described by the 4-vector poten- 
tial = (Aq,A) with components related to its rest 
Coulomb potential U(r): 



A =-yU A = j/3U 



(1) 



Here /3 = v/c, v is the target atom velocity in the EA 
rest frame and 7 is its Lorentz-factor. The time-like com- 
ponent Aq of the 4-potential interacts with the charges of 
the particles forming EA and the space component with 
their currents. 

In this paper we consider only EA consisting of spinless 
particles (tt, if-mesons etc.) which are of interest for 
the DIRAC experiment. In the Born approximation the 
amplitudes of transition from the initial state i to the 
final one / due to the interaction with A^ can be written 
as: 



Afi = U{Q)a fl {q) , 

OO 

U(Q) = 2 J U(r)^rdr, 


a Mi) = pfi{i)-Pjfi(Q) ■ 



(2) 
(3) 
(4) 



The transition densities Pfi(q) and transition currents 
jfi(q) are expressed via the the EA wave functions tpi 
and tj)f for the the initial and final states: 



" „—iq-ir \ j3 



Pfi(Q) = J P fl {r)(e^ r -e~^ r )d?r, (5) 



d A r , (6) 
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Pfi (r)=r f (rWi(r) (7) 
[Mr)Vr f (r)~r f (r)VMr)] ■ (8) 



The EA wave functions ipij and the binding energies 
£ij obey the Schrodinger equation 



Hipij = Sijipij 



(9) 



with the Hamiltonian H. It is worth noting that the 
explicit form of the potential V(r) of the interaction be- 
tween the EA components have no influence on the final 
result of this paper. 

In the above equations m 12 are masses of EA com- 
ponents, q — (qo,g) is the transfer 4-momentum. All 
other kinematic variables are related by the following 
equations: 



<7i 



JL 
mi 



A* 

= — Q , M 

m 2 



= M 

q = (qL,qT 



M = mi + 7712 



q = uj fl + — = (3q = (5q L . 



(10) 



Q = VQ 2 



Q 



•d(i 



P 2 



The differential and integral cross sections of the EA 
transition from the initial state i to the final state / due 
to interaction with the electromagnetic field of the target 
atom are related to amplitudes (J2j): 



CTfi 



1 

7p 



\A f M\ 2 d 2 q T ■ 



(11) 



Formulae |TT|) allow to calculate the transition (par- 
tial) cross sections in the Born approximation. But for 
applications (for example see Q) the total cross sections 
of the EA interaction with target atoms are also required. 
Because the Born amplitudes of the EA elastic scattering 
are pure real values, the optical theorem cannot be used 
to calculate the total cross sections. Thus they should be 
calculated as the sum of all partial cross sections: 



tot \ " 

f 



(12) 



To get a closed expression for the sum of this infinite se- 
ries (the so-called "sum rule" ) the transition amplitudes 
W) are usually rewriten as: 



A fi (q) = (f\A(q)\i) , 



(13) 



where the operator A(q) does not contain explicit depen- 
dence on the EA final state variables (for example, its 



energy Ef , see bellow). Then, using the completeness 
relation 



E i/></i = 1 > 

we can writte sum dl2) in the form: 



(14) 



/3 2 



(i\A*(q)A(q)\i)d 2 q T . (15) 



III. SIMPLIFIED APPROACH 

One should take some caution when passing from the 
exact expressions (p|~[lO|) for the transition amplitudes, 
with explicit dependence on the e f (through the time- like 
qo and longitudinal q^ components of 4- vector q), to the 
approximate one without such dependence. Otherwise, 
it is possible to obtain a physically improper result as it 
has happened to the authors of the paper [Q at deriving 
of the sum rules for the total cross section of interaction 
of ultrarelativistic EA [j3 = 1) with target atoms. Below 
we discuss this problem in detail. 

The most essential simplification, that arises in the 
case of = 1 is that Q 2 = q\. Thus U(Q) = U(q T ) 
[see (10)] and only Afi in (|2|) depends on Sf through the 
exponential factors exp (iq±r) and exp {—iq%r) in (m and 

(I) 



<?1,2T- = 



mi,2 



-qr 



"11.2 



(qLZ + qrrr) 



(16) 



where qh — w/j + q^/2M if /3 = 1. 

Now let us take into account the fact that the typical 
value of z in these expressions is of the order of the Bohr 
radius tb — 1/fia and the typical q^ ~ ujf i ~ l^a 2 , thus 
the product q^z is of the order of a. Then it seems 
natural to neglect the g^-dependence of ajf. 



a fi (q) w a fi (q T ) 



(17) 



and consider this case as the zero order approximation 
to the problem It corresponds to the choice of the 
operator A in the form: 



A(q) = U(q T ) [e*» 

(e^ TrT /mi + e - iq2TrT /m 2 )/3p] 

Here p = — iV is the momentum operator. 



(18) 



Substituting (18) in ( |15| ) results in the following sum 
rules H, where the total cross section is expressed as 



the sum of the "electric" cr cl and "magnetic" <r r 
sections: 

tot el I _mag 

it = a + a . 



a cl = J U 2 (q T )M(q T )d 2 q T 

M(q T ) = 2(1 - S(q T )) 
S(q T )= f mr)\ 2 e^ r d 3 r: 



cross 

(19) 
(20) 



3 



K(q T ) = J 



J U 2 {q T )K{q T )d 2 q T , 



(21) 



1 



/i z 77117712 



-{e lqr - 1) 



|/3 V Mr)\ 2 d 3 r 



These results reproduce the ones obtained in m and 
differ from the sum rules used in ^ by the additional 
term <7 mag . For beginning let us consider its contribution 
qualitatively. For this purpose the target atom poten- 
tial U(r) can be approximated by the screened Coulomb 
potential: 



TTI \ Za -Ar 
U (r) = e 

r 



A 



(22) 



where m e is the electron mass and Z is the atomic num- 
ber of the target. The pure Coulomb wave function can 
be used for ipi (i.e. the contribution of the strong interac- 
tion between the EA components is neglected, see @). 
For the ground state it is written as: 



no?!" 2 



(23) 



Under such assumptions for the ground state the follow- 
ing results can be easily obtained: 



mag 



8vrZ 2 , 

5— m 

M . 
4tt (Za 
T \T 
ATtZ^a 2 
3m 2 



2fi 



Z^l^rrie 



+ 0(a 2 a cl ) 
+ 0{a 2 a cl ) . 



(24) 



(25) 



It is seen that in spite of a 2 in the numerator of <r mas 
the electron mass square in the denominator makes the 
contribution of the "magnetic" term in (|l^) not negligible 
with respect to the "electric" one, especially for the case 
of EA consisting of heavy hadrons and low Z values. 

To obtain exact numerical values we have precisely 
repeated the calculations made in jTj. More accurate 
presentation of the target atom potential, namely, the 
Moliere parametrization of the Thomas-Fermi potential 
p8[ was used as in : 



U(r) = ZaJ2 



ae- x ' r 



(26) 



ci=0.35, c 2 = 0.55, c 3 = 0.1; 
Ai = 0.3A , A 2 = 1.2A , A 3 = 6A , A = m e aZ 1/3 /0.885. 

The values of the "electric" (el) and "magnetic" (mag) 
total cross sections (in units of cm 2 ) and their ratio 
(mag/el) are presented in Table [j| for various EA and 
target materials. The values published in Q are given in 
parentheses. It is seen that the "electric" cross sections 
coincide within the given accuracy, but the "magnetic" 
ones are underestimated in Q. It is worth noting that 



the correct values of c mag do not depend on EA masses as 
it follows from the simplified approximation result j2^). 
The ratio values confirm the above estimation about the 
"magnetic" term contribution. Thus, inaccuracy in the 
calculations did not allow the authors of JtJ to observe 
such a significant contribution of cr mas in their results. 



TABLE I: The "electric" (el) and "magnetic" (mag) total 
cross sections in units of cm 2 and their ratio (mag/el) in % for 
EA consisting of tv and K mesons {A^-n, A^k, A^k) and target 
materials with the atomic number Z. The values published 
in |7f are given in parentheses. 
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It is clear that such strong enhancement of the mag- 
netic term in (Il9[) is the consequence of its inverse power 
dependence (|25f) on the small screening parameter A. It is 
also easy to see that the origin of such unnatural depen- 
dence is in the behaviour of the factor K(qx) at small qr 
in (pl|). This factor, contrary to M(qt) in (p0|), does not 
approach zero at q? — * 0. But at (3 = 1 such behaviour 
of K(qr) is in conflict with some general properties of 
transition amplitudes (Q), which follow from the continu- 
ity equation: 



Ufipfi(q) ~ qjfi(q) = 



(27) 



(The latter can be derived from the Schrodinger equation 
(0)). Indeed, rewriting the continuity equation in the 



4 



form: 

wfiPfi(q) - qLPjfi(q) - qrjfiiq) = 

Ufiipfi ~ PJfi(q)] - qlPjfriq)/™ - q T jfi(q) = , 

(28) 

it is easy obtain 

a fi(q) = pfi(q)~ Pjfi(q) 

1 (29) 
= —[qlP3fi{q)/1M + q T j fi {q)]. V > 

Thus all transition amplitudes become zero at qx = 0. 
Therefore, any transition cross section ( pi] ) can depend 
on the screening parameter A at least only logarithmi- 
cally, but never like inverse power of this parameter. The 
same is valid for the sum ( |l2|) of this quantities, i.e the 
total cross section. 



IV. ACCURATE FORMULAS 

Since the A-dependence of the magnetic term in fl25| ) 
is contradictory to the general result, we must conclude 
that there is a fallacy in the deriving of sum rules @ 
somewhere. To understand the origin of the error, made 
by the authors of J7|], let us go back to quantities (||),@ 
and expand them in powers of the longitudinal momen- 
tum transfer q^'- 



Pfi 



Pf 

71 = 



(n) 



in) _Ql( d n 



E.(«) 
3 }i ■ 

n=0 



,•(«) 
J ft 



<fl ( dn 



n\ \dql Jfi 



q L =0 



9i=0 



(30) 



(31) 



It is easily shown that terms of these expansions obey 
the following estimation: 



(") r, 

P)i' cx a 



(32) 



The additional power of a in the current expansion coef- 
ficients, in comparison with the density one, reflects the 
ordinary relation between the values of current and den- 
sity in the hydrogen-like atoms. 

Expanding (|J) and taking into account ([52]) it seems 
reasonable to group terms with the same order of a rather 
than q L as was done in [fjj. Then the successive terms of 
the afi expansion in powers of a are 



a fi 



p)7 & 



(33) 



fi 



From above it is clear that in the "natural" approxi- 
mation ( |l7| ) includes a$ and only one part of the term 



"fi 

of expansion (p3|), namely: 



(34) 



while the second one 



pfi=iqL J 4>fE(q T ,r T )zipid 3 



(35) 



was omitted according to the reasoning of approximation 
(|17|). In equations (BJ), (j35|) E(qT,rx) denotes: 



E{q T ,r T ) 



P t 



P 
( 

m 2 



(36) 



Let us consider this neglected part in detail. As it is 
proportional to qL = oJfi + q^/2M and therefore explic- 
itly depends on Ef, one cannot use completeness relation 
( |l4| ) to calculate its contribution to the total cross sec- 
tion directly. First we need to transform it to the form 
free of such dependence. It can be done with the help of 
Schrodinger equation (||). 



s f i J ip* f {r)E{q T ,r T )zip l {r)d i r = 
J ip)(r) {e f E(q T ,r T )z - e l E{q Tl r T )z}^ l {r)d 3 r = 

r f (r)[H,E(q T ,r T )zWr)d 3 r (37) 



The commutator in this relation is easily calculated and 
after simple algebra we get the following result: 



i>}(r)E(q T ,r T ) 



dipj(r) 
dz 



(38) 



Apg(<?) = 



mi 



rn-2 



Ox, 



q\ ± 2q T p 

2mi.2 



zipi(r)d 3 r 
(39) 

iV . (40) 



It is seen that "large" (nonvanishing at qr = 0) parts 



of two terms (34) and (|38j), contributing to , are equal 
and opposite in sign, so that in the resulting expression 
they cancel each other, leaving only the term with the 
"correct" behaviour at small q^: 



x (1) 



» 



(41) 



The same is valid for any . Applying Schrodinger 
equation @ to reduce by one the power of q^ in the 
expression: 



p f 7(q) 



(iQLY 



Mr) 



+ (-1) 



n+l 



m 2 



P 

mi 

-iq2TfT 



a iqiTTT 



z n Mr)d 3 r , (42) 



one can represent it in the form: 



pfhq)^(3j ( ;r 1 \q) + ^{q) (43) 



5 



A^)( 9 )=^ 



mi 



m 2 



z n if>i(r)d 3 r 



Finally one has 



An) 



(44) 



(45) 



That confirms the qualitative result (29), derived with 
help of continuity equation (27). 



The remaining e^-dependence of right-hand side of 
(44|) can be removed by repeatedly applying the Schro- 
dinger equation (^|), which allows the transition ampli- 
tudes to be represented in the form (ft3|). 

From ^-dependence of the integrand in (|44|) it is easy 

to derive that afp = for odd A(lm) jj, and a^* = 
for even A(Zm)/j , where A(Zm)/j — (If — li) — (my — 
mj), and li, If, nii, nif are the orbital and magnetic 
quantum numbers of the initial i and final / states (the 
quantization axis is supposed to be z-axis). Thus "odd" 
and "even" terms of expansion ( |33| ) do not interfere and 
therefore in the expansion of the cr tot in the powers of a 



oo 
n=0 



a {n) cx a % 



(46) 



only even powers are present. 

The structure of the zero order term of this expansion 
is well established [see (pp|)]. In view of the above dis- 
cussion one can be assured that the higher order terms 
are numerically negligible and should not be discussed 
in detail. Nevertheless, for completeness of the consid- 
eration we present the expression for the contribution of 
the a 2 -term to the total cross section which includes the 



term 



l fi 



and the interference term a ?^ . 



r (2) 



U 2 {q T )W(q T )d 2 q T + 0(a 4 ) 



(47) 



W(q T ) = 



1 



4mim 2 



z 2 [q"mrt 



\2q T pipi{r)\' 



lqr d?r 



The "correct" gr-dependence of the last integrand ex- 
cludes a possibility of some extra A-dependence arising, 
which could dramatically enhance the contribution of this 
term (as happened to the cr mas term in jjj). This can 
be illustrated by the explicit expression for the case of 
the screened Coulomb potential (|2^) and the EA ground 
state (|23|): 



a< 2 ) = - 



&n{ZaY 
5Mfi 



In 



2//. 



ZV3, 



(48) 



Because of numerical smallness of a 2 this term can be 



successfully neglected compared to (J24J) . 

Thus in most practical applications, in which the re- 
quired relative accuracy is less than 10 -4 , only the zero 
order term, that considers the pure Coulomb interaction 
and only the transverse transfer momentum, should be 
taken into account for calculation of the relativistic atom- 
atom cross sections. This result warrants the usage of the 
simple expression: 



U 2 {q T )[l~S(q T )}d 2 q T 



(49) 



for the total cross section calculation for the Born ap- 
proximation in Q and for the Glauber extensions in Jl5[ . 

The authors express their gratitude to professors 
S.Mrowczyhski, L.Nemenov and D.Trautmann for help- 
ful discussions. This work is partially supported by 
RFBR grant 01-02-17756. 



[1] B. Adeva et al., Lifetime measurement of-K + n atoms to 

test low energy QCD predictions, Proposal to the SPSLC, 

CERN/SPSLC 95-1, SPSLC/P 284, Geneva, 1995. 
[2] L.G.Afanasyev and A.V.Tarasov: Yad.Fiz. 59, 2212 

(1996); Phys.Atom.Nucl. 59, 2130 (1996). 
[3] G.H.Gillespie, Phys.Rev. A18, 1967 (1978); 

G.H.Gillespie, M.Inokuti, Phys.Rev. 22A, 2430 (1980). 
[4] L.S.Dul'yan and A.M.Kotsinyan: Yad.Fiz. 37, 137 

(1983); Sov.J.Nucl.Phys. 37, 78 (1983). 
[5] A.S. Pak, A.V. Tarasov: JINR-P2-85-903, Dubna 1985. 
[6] S.Mrowczyriski: Phys.Rev. A33, 1549 (1986). 
[7] S.Mrowczyriski: Phys.Rev. D36, 1520 (1987); 

K.G.Denisenko and S.Mrowczyriski: Phys.Rev. D36, 

1529 (1987). 

[8] A.V.Kuptsov, A.S. Pak and S.B.Saakian: Yad.Fiz. 50, 

936 (1989); Sov.J.Nucl.Phys. 50, 583 (1989). 
[9] L.G.Afanasyev: JINR-E2-91-578, Dubna 1991. 



[10] L.G.Afanasyev and A.V.Tarasov: JINR E4-93-293, 
Dubna, 1993. 

[11] L.G.Afanasyev: Atomic Data and Nuclear Data Tables, 

61, 31 (1995). 
[12] Z.Halabuka et al.: Nucl.Phys. 554, 86 (1999). 
[13] A.V.Tarasov and I.U.Christova, JINR P2-91-10, Dubna, 

1991. 

[14] O.O.Voskresenskaya, S.R.Gevorkyan and A.V.Tarasov: 
Yad.Fiz. 61, 1628 (1998); Phys.Atom.Nucl. 61, 1517 

(1998) . 

[15] L.Afanasyev A. Tarasov and O.Voskresenskaya: J.Phys. 

G 25, B7 (1999). 
[16] D.Yu.Ivanov and L.Szymanowski: Eur.Phys. J.A5, 117 

(1999) . 

[17] Amirkhanov I. et al., Phys. Lett. B 452, 155 (1999). 
[18] Molfere G., Z. Naturforsch. 2A, 3 (1947). 



